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Abstract 



In this paper we present a detailed formulation for a recently proposed ef- 
fective field theory to describe the nonperturbative QCD dynamics of heavy 
mesons. This effective theory incorporates with heavy quark symmetry (HQS) 
and the heavy quark effective theory (HQET). Heavy mesons in this theory 
are constructed as composite particles of a heavy quark bounded with the 
light degrees of freedom. The heavy meson properties in the heavy quark 
limit and the I/tuq corrections can then be explicitly evaluated from this 
effective theory. All the basic parameters of the HQET, namely, the heavy 
quark mass uiq, the heavy meson residual mass A, and the HQS breaking 
mass parameters Ai and A2, are consistently determined. Ai is found to be 
small due to a large cancellation between the heavy quark kinetic energy and 
the chromo-electric interaction between the heavy quark and light degrees of 
freedom. We also evaluate the Isgur-Wise function, the decay constant, and 
the axial-vector coupling constant of heavy mesons. 



PACS numbers: 12.38.-t, 12.39.Hg, 12.39.Ki, 12.60.Rc, 14.40.-n 
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I. INTRODUCTION 



Within the last decade, the most significant progress in the QCD description of hadronic 
physics is the discovery of heavy quark symmetry (HQS) which, together with the de- 
velopment of heavy quark effective theory (HQET) from QCD ||, has largely simplified the 
analysis of heavy hadron physics ||. However, even in the infinite quark mass limit, the 
general properties of heavy hadrons, namely, their decay constants, transition form factors 
and structure functions etc. measured in various exclusive and inclusive decay processes re- 
main unknown in the framework of QCD. The difficulty in understanding the heavy hadronic 
properties essentially lie in nonperturbative QCD dynamics, but HQET itself does not pro- 
vide such a description to nonperturbative QCD. Lattice QCD simulations 0] permit a 
nonperturbative approach to the low-energy QCD problem, but so far a direct lattice calcu- 
lation with heavy quarks is still not possible due to the difficulty of placing heavy particles 
on the lattice. An alternative first-principles calculation of nonperturbative QCD to these 
heavy hadron properties is to solve the heavy meson bound states directly from light-front 
QCD ||, but to obtain reliable results, further investigations are needed. Therefore, in most 
recent studies, these heavy hadron properties are usually evaluated using phenomenological 
models, such as the constituent quark model ||, the MIT bag model 0, the QCD sum rule 
||, and the light-front quark model 0. 

As is well known, although the constituent quark model (CQM) and the MIT bag model 
have been widely used in the phenomenological discussion of hadronic structures, appli- 
cations of these two models are trustworthy only for processes involving small momentum 
transfers. This is due mainly to the nonrelativistic limitation of the CQM and to the difficulty 
with boost in the bag model. The light-front quark model (LFQM) which is a relativistic 
quark model with simple boost kinematics allows us to describe physical processes with 
large momentum transfers, but it is still not truely Lorentz covariant due to the exclusion of 



the so-called Z-diagrams |T0|,|TTJ. As a result, certain theoretical ambiguities arise in LFQM 
calculation which may lead to inconsistent results as have been shown in [p~2| , |13|| ■ 

To overcome the drawbacks in these phenomenological models, we recently proposed a 
covariant light-front model JH] which modifies the conventional LFQM in the heavy quark 
limit by adding a constraint on the light-front wave function. The covariant light-front 
model rules out some non-covariant light-front wave function often used in the literature, and 
therefore partially eliminates the ambiguities presented in previous calculations. Meanwhile, 
it also largely simplifies the light-front formulation, and may further provide a first-principles 
QCD analysis of the 1/rriQ corrections within HQET ||T4 |. 

However, further investigation indicates that although it has overcome some theoretical 
difficulties encountered in the conventional light-front formulation for heavy hadrons, our 
covariant quark model |14| still cannot provide fully consistent results for processes involv- 
ing light quark currents. This is because for those processes involving only heavy quark 
transitions, the Z-diagram contributions are suppressed in the heavy quark limit; however, 
when the light quarks (or currents) are involved, the non-covariant light-front treatment of 
light quarks (due to lack of Z-diagram contributions from the light quark production) will 
still cause theoretical ambiguities. 

In this paper, we will present a detailed formulation of a fully field theoretical description 
of heavy mesons in terms of an effective theory of composite particles we have proposed very 
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recently |L5| . The composite particle consists of a reduced heavy quark (heavy quark in the 
heavy quark limit) coupled with the light degrees of freedom, in which a structure function of 
the composite particle, ty(v-p q ), corresponding to the wave function of a heavy meson bound 
state, is explicitly built in. In this field theory description, the constituents in composite 
particles are no longer fixed at a given instant point of time, which is a condition imposed 
in the usual construction of hadronic bound states in various quark models. The field 
theory structure of the constituents in hadrons allows us to formulate the physical processes 
in terms of the standard Feynman diagrams in which various time-ordering diagrams are 
all automatically included. Therefore, the lack of Z-diagrams in the usual quark model 
descriptions is no longer a problem in the present formulation. 

Moreover, combining the effective Lagrangian of the composite heavy mesons with the 
expansion of the heavy quark QCD Lagrangian, we can systematically evaluate various 
1/m.Q corrections to heavy meson properties in the standard framework of perturbative 
field theory. Thus, a self-contained description of heavy mesons (including the bound state 
structure and I/tuq corrections) is realized in a field-theoretic framework. This effective 
field theory allows us to explore the nonperturbative heavy meson dynamics, which is not 
possible in the light-front quark models. 

The rest of the paper is organized as follows: In Sec. II, we will analyze the basic structure 
of heavy mesons in the heavy quark limit as a composite particle of the reduced heavy quark 
coupled with the light degrees of freedom. In Sec. Ill, we construct an effective Lagrangian 
to describe the composite structure of heavy mesons and combine this effective Lagrangian 
with HQET to establish a realistic effective field theory which can be used to evaluate the 
1/rriQ corrections in the standard Feynman diagrammatic approach. As first applications, 
we evaluate the Isgur-Wise function, the heavy meson decay constant and the axial-vector 
coupling constant in the heavy quark limit in Sec. IV. In Sec. V, we compute the pseudoscalar 
and vector heavy meson masses up to 1/toq, which determine the HQET parameters Ai and 
A2. In Sec. VI, we compare the effective field theory with the covariant light-front model 



14| , and show how the effective field theory overcomes the lack of relativistic covariance in 
light-front quark models. In Sec. VII, we present some numerical calculations to check the 
self-consistency of the theory, and then calculate all the basic parameters in HQET. The 
summary and perspective are given in Sec. VIII. 



II. A COMPOSITE PARTICLE PICTURE OF HEAVY MESONS IN THE HEAVY 

QUARK LIMIT 

We define the general expression for the composite operators of pseudoscalar and vector 
heavy mesons as follows: 

H ci (X) = J d 4 yq(X - ay)Y i Q[X + (1 - a)y] , i = P, V , (2.1) 

where the subscript c means "composite", Q and q are the heavy and light quark field 
operators respectively, X is the center-of-mass coordinate of the heavy meson and y = xq—x q 
the relative coordinate between the heavy and light quarks, a = itiq/ '(mQ+m q ), and Tp = 75 
and Ty = 7 M define the spin structures for the pseudoscalar (0~) meson and vector (1~) 
meson, respectively. Inside the heavy meson, QCD dynamics is nonperturbative, so that 
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Q(x) and q(x) are strongly coupled, and they are surrounded by infinite number of qq 
pairs and gluons originating from the nontrivial QCD vacuum. We may phenomenologically 
rewrite the above composite operator of heavy mesons in terms of the constituent valence 
heavy and light quark field operators, denoted by Q and q respectively, coupled through 
a structure function F(y) which describes the binding effect of infinite number of qq pairs 
and gluons governed by QCD, 

H Ci (X) = J d 4 yq (X - ay)F(y)T i Q [X + (1 - a)y] . (2.2) 

Now we take the heavy quark limit, i.e., rriQ — > oo. In the momentum space, the heavy 
meson carries momentum P M = rriQV M + p H , where t> M is a four velocity (v 2 = 1) of the 
heavy meson, and pn the residual momentum. The heavy meson on-mass-shell condition 
P H = Mhv^ corresponds to p H = Av 11 , where Mh is the heavy meson mass which approaches 
to infinity under the heavy quark limit but A = Mh — rriQ is kept finite. Note that A is a 
basic parameter in HQET known as the residual mass of the heavy meson. 

In the heavy quark limit, we define the so-called reduced heavy quark field h v in the 
heavy quark expansion ||, 

Q (x)=e- im ^ x h v + O(l/m Q ) , ±±lh v = h v . (2.3) 
We also introduce the reduced heavy meson composite operator H ci , 

H ci (X) = -±^e- im ^ x H ci (X) . (2.4) 

V M H 

Then the Fourier transformation of the composite field operator in the momentum space 
can be expressed as 



H«(v,p H ) = Jd*Xe^ x H a (X) 



(2tt) 4 (2tt) 4 ' 

xqM^iv-p^TMk), (2.5) 

where = Pq — vtiqv^ is the residual momentum of the heavy quark, and p q the momentum 
of the light antiquark. The function ^>(v ■ p q ) is the heavy-quark-limit expression of the 
Fourier transformation of F(y), which is the analog of heavy meson wave function in the 
covariant light-front quark model |14| . 



One may ask why is ty(v-p q ) only a function of v-p q ? Note that the Fourier transformation 
of F(y) should be a scalar function of the relative momentum q = Pq — p q = rriQV + k — p q , 

F(q 2 ) = ( d A ye^F{y) . (2.6) 



In the heavy quark limit, the k - and ^-dependences are suppressed by the heavy quark 
mass: 

q 2 = m Q (m Q + 2v(k- p q ) + 0(l/m Q )) . (2.7) 
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Also, due to heavy quark symmetry, F(q 2 ) must be a m^-independent function in the heavy 
quark limit. Besides, by momentum conservation [i.e. the 5-function in Eq.( |2.5| )], k ~ Pu—Pq 
and hence q 2 ~ (v ■ pu — 2w • p q ) . Since \1/ describes the composite structure of heavy meson 
bound states in which pn = Av, then, without loss of the generality, we have replaced F(q 2 ) 
by the mQ-independent function \l/(t> ■ p q ) in Eq. ([2.5|). 

A motivation of the above analysis is to provide a realistic field theory formulation of 
the intuitive picture of heavy hadrons; that is, in the heavy quark limit the heavy meson is 
a composite particle which can be described by a composite field operator of a heavy quark 
h v (k) coupled with the light degrees of freedom which is sometimes called the "brown muck" 
in the literature. From Eq. ( |2.5| ), we see that it is natural to define the brown muck by 

q v {p q ) = V{vp q )q {p q ). (2.8) 

The above definition indicates that the brown muck consists of a light valence antiquark 
(which contains the tensoral structure of a spin-1/2 Dirac particle) coupled with a brown 
muck structure function ^(v-p q ), which effectively describes the dynamics of infinite number 
of qq pairs and gluons surrounding the light antiquark. Thus, according to Haag, Nishijima 
and Zimmermann |T7j], we may rewrite the composite field operator of heavy mesons as a 
local operator: 

H ci (x) = q v (x)Tih v (x) , (2.9) 
and its momentum representation is given by: 

H a (v,Pn) = J (0i^( 27r ) 4 ^(^ - k - P q )*MTMk) . (2.10) 

From Eq. ( j2.4| ), the normalization condition of the heavy meson bound states (with pn = 
Av ) in the heavy quark limit is given by 

(H ci ,(v')\H ci (v)} = 2(2vr)V<5 3 (Ai; - Av')5 u , , (2.11) 

where the superscript e denotes the energy component in the momentum space, which can 
be either the or the + component, depending on whether the light quark on-shall energy 
is picked on the light-front, p~ = (p 2 x +mjj)/p+, or on the equal-time form, p° = ^Jp^ + m 2 . 

Correspondingly, 5 3 (Av — Av') becomes 8{Av + — Av' + )5 2 (Av±_ — Av' ± ) or 5 3 (Av — Av'). 

Of course, a first-principles determination of the above composite particle picture lies in 
the detailed form of \l/(f ■ p q ). How to solve ty(v ■ p q ) directly from QCD is one of the most 
interesting and difficult problems in strong interaction physics. In the present paper, we 
shall treat \l/(f ■ p q ) in a phenomenological manner. 

III. EFFECTIVE FIELD THEORY OF HEAVY MESONS 

Based on the above analysis of the composite particle picture of heavy mesons in the 
heavy quark limit, we can now build an effective field theory to describe the heavy meson 
structure in HQET. 
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A. Effective Lagrangian 



In principle, the heavy-meson composite particle structure is determined by the QCD 
Lagrangian for heavy and light quarks plus gluons, 

C = Cq + C q + Cg 

= QiilP - m Q )Q + q{iip - m q )q - -F^ 'F a/Xu . (3.1) 

Using the I/tuq expansion to the heavy quark QCD Lagrangian ||, the above Lagrangian 
can be rewritten as 

C = Cq + C mQ , (3.2) 

where 

C = h v iv ■ Dh v + q(ip - m q )q - ^F^F aiMU (3.3) 



governs the structure of the heavy meson composite field operator of Eq. (|2.9| ), and it has the 
heavy quark ST// (2) ® SU S (2) flavor-spin symmetry (or simply the heavy quark symmetry, 
HQS), while, 

C m Q = E ( — J h v (ifl±)(-iV ■ D) n -\iP ± )h v 
oo / i \ n 

(3 - 4) 

determines the I/tuq corrections to Cq, which breaks the heavy quark symmetry, where 
D^_ = D>> - v^v ■ D. 

To directly solve the structure of the composite field operator of Eq. (|2.9| ) from C is not 
simple and may not even be possible within the known framework of field theory. Instead 
we shall introduce an effective Lagrangian to phenomenologically describe the heavy meson 
composite structure, 

C -> C^(h v , q, H ci ) = h v iv ■ dh v + q (i - m q )q 

+Pl(iv d -2A )P V - Vf(iv d -2A )^ 
+G (h v i'j 5 q v P v - HvjpqvVf + /i.e.) , (3.5) 

which has the same heavy quark symmetry as Co, where d=d — d, q v is the brown muck field 
of the light degrees of freedom inside the heavy mesons, specified by Eq. ( |2.8[ ) in momentum 
space. P v and V£ are the reduced pseudoscalar and vector heavy meson fields, respectively: 

$(x) = -^e-^^(x), (3.6) 

A*(x) = -±=e*»<>**Vf(x), (3.7) 

V My 
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with v ■ V v — 0, and Ao = Mq — rriQ as the bare residual mass of the heavy mesons, where Mq 
and m,q are defined here as the bare masses of heavy mesons and heavy quarks, respectively. 
The third and fourth terms in ( |3.5|) are directly obtained from the free Lagrangian of the 
pseudoscalar and vector fields with the above definition of the reduced fields, 

C ee = - M|,|$| 2 ) - (^F^F, U + M^A^A,) , (3.8) 

where F» v = d»A u - d v A». 

In principle, the effective Lagrangian may be derived from Cq by integrating out the 
gluon degrees of freedom in a nonperturbative way, which is unfortunately not practical at 
present. Therefore, we shall assume that after integrating out the gluon degrees of freedom 
mediated between heavy and light quarks, we are led to an effective Lagrangian of the form: 

£gff = h v iv ■ dh v + q o (i0- m q )q Q 

+go\h v ij 5 q v q v ij 5 h v - hvl^q^K H ) , (3.9) 

in which the complicated nonperturbative QCD dynamics to the heavy meson structure is 
effectively described by introducing the brown muck field q v through the universal structure 
function \l/(f • p q ) of Eq.( |2.8| ). A consistency condition for the composite particle picture of 
heavy mesons in the heavy quark limit can be determined by demanding the equivalence of 
Eqs. (\3.d[) and ( \3.£\ ). Note that a priori the structure of q v (x) can be different in Eqs.(|3.5|) 



and (|3.9| ), however the consistency requirement demands that they must be the same fll5| . 



B. Dynamical description of the composite particle structure 



We can now dynamically describe the heavy mesons in the heavy quark limit as composite 
particles in the above two effective field theories. In terms of Eq. ( |3.9|) , the composite particle 
structure of a heavy-light quark field q v Yih v can be obtained by considering the heavy-light 



quark scattering in the chain approximation, as shown in Fig. 1 (T^|. In the pseudoscalar 
channel, apart from the factors come from the external quark lines, the scattering amplitude 
is given by 



(3.10) 



where pn = k +p q is the momentum transfer, k the residual momentum of the heavy quark, 
p q the momentum carried by the light quark, and n#(f • pu) the "self-energy" correction 
from the heavy-light quark loop depicted in Fig. 1(a), 



- l U H (vp H ) = (^r(-l) J ^Mv-p q )\ 2 Tr 



(27T) 



T P i 



1+jf 



2(v ■ p H - v ■ p q + ie) 



x Tpi 



- i> q + m q 
v\ - m \ + * e 



(3.11) 



where Tp = 275. Thus, we have 
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TT (v n \ - i [ ^ \^( V -Pi)\ 2 2{vp q + m q ) 

LL H {V -p H ) — if — -— — — — - , [6.LZ) 

J {2irp (v ■ p H - v ■ p q + le) (p* - + ie) 

The existence of a stable heavy pseudoscalar composite meson with the residual mass A 
implies that the above amplitude has a pole at pn = (which corresponds to Ph = Mhv). 
This leads to 

= (313) 

If we expand Hh(v • Ph) by 

K H {v ■ p H ) = riff (A) + {vp H - A)n' H (A) + U r H {v ■ p H ) , (3.14) 
then, the amplitude can be expressed as 

, y(vp q )V*(vp' q ) 

Qq ~ \v.p H -A)W H (A) + U H (vp H )- [6Ab) 



On the other hand, in Eq. ( |3.5| ) the heavy meson field appears as a fundamental field. 
The physical (i.e., renormalized) meson structure can be determined by considering a similar 
heavy- light quark scattering process shown in Fig. 1(b). The amplitude is 

Am = Gfo{v ■ p q W{v ■ p' q )A H (v ■ p H ) , (3.16) 

where A(v ■ pn) is the heavy meson propagator: 

A h(v ■ Ph) = -, =f r * 2n , r , (3.17) 

2(v ■ p H - A ) - G$l H {y ■ p H ) 

with Tl H (v ■ pn) being given by Eq. ( |3.12| ). Using the expansion (|3.14|) , the above amplitude 
can be recast into 

AM = G 2 *(v Pq )**(vp' q )- - — -, (3.18) 

2(v ■ p H - A) - G 2 W H {v ■ p H ) 

with 

G = Z 1 3 /2 G , Z 3 = l + ^-U' H (K), (3.19) 

A = A + ^W(A), (3.20) 

where Z3 is the wave function renormalization constant of the heavy meson field P v . 

In order that the pseudoscalar heavy meson, being a composite particle with the struc- 
ture q v j5h v , can legitimately be represented by the field operator P v , the two scattering 
amplitudes Eqs. ( |3.18| ) and ( |3.15| ) must be the same. This results in 

,-2 . l TT > nr\ ■ f \^( v -Pq)\ 2 v-p q + m q 



G - 2 U' H (A)-ij {27r ^ (K _ v . pq 9 +K)2 {p2 _ q m 2 + l e y ( 3 - 21 ) 
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and 



A = A + G 2 U H (A) 

_ x i ^ f d A p q \^(v ■ p q )\ 2 vp q + m q 

- A + ZC J ( ^ )4 ( _ _ v ^ + ^ _ ^ + . £) • l^J 

We therefore obtain from Eq. (|3.19|) that 

Z 3 = . (3.23) 

The above results have a clear physical interpretation. The fact Z3 = implies that the 
bare fundamental field P v = Z\J 2 P R does not exist. In other words, the physical meson 
field P R must be a composite particle operator. Such a realization of the composite particle 
relies upon the introduction of the composite particle structure function ty(v ■ p q ) in the 
effective Lagrangians (|3.5|) and ( |3.9|) . One can always choose a structure function ty(y ■ p q ) 
to ensure the finiteness of II#(A) and IT^ (A). Without invoking such a structure function 
in the effective four-fermion point interaction, T\.' H (Qr) is divergent and hence G = 0, which 
leads to obvious inconsistency in the effective theory. This shows the importance of ^>(v-p q ) 
in our construction, which is not surprising since, as we will see shortly, $?(v ■ p q ) actually 
corresponds to a heavy meson wave function, while the renormalized coupling constant G 
of Eq. (|3.21|) is just the wave function normalization constant. 

Similar discussion for the vector meson structure can be easily carried out. In the heavy 
quark limit, the heavy quark symmetry of the effective theory ensures that the vector meson 
composite particle has the same structure as the pseudoscalar particle. We will therefore 
not repeat the similar derivation here for vector mesons. 

C. The effective field theory of heavy mesons with Feynman rules 

After determining the composite particle structure of heavy mesons in the heavy quark 
limit, we can proceed to evaluate various heavy meson properties by combining the effective 
Lagrangian C^ R , which obeys heavy quark symmetry, with the 1/rriQ corrections from C mQ : 

£cS = £cS R + An Q , (3.24) 

where 

C% R = h v iv ■ dh v + q(i$- m q )q 



+Pl{iv d -2A)P V - V v rt(iv d -2A)V Vfl 

+G(h v ij 5 q v P v - Ki^V* + h.c) , (3.25) 
00/ 1 \ n _ 

£m Q = E J h v {ifl ± ){-iv ■ D) n ~ x (i]p}_)h v , (3.26) 

G is determined by Eq. ( |3.21| ), A is the physical residual mass of heavy mesons in the 
heavy quark limit, M H = tuq + A, q v = "$?(v ■ p q )q and ^f(v ■ p q ) describes the heavy meson 
structure and is a phenomenological input at this level. The Lagrangian C^ R gives the 
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nonperturbative effective coupling of heavy mesons with heavy-light quarks. And C mq is 
treated as a perturbation to C^ R , which contains all the I/ttlq corrections. Then a practical 
evaluation scheme can be developed in terms of the standard Feynman diagrammatic rules: 
(i) The heavy meson bound state (p# = Av) in the heavy quark limit gives a vertex as 
follows: 

<^ : -iGV*(v Pq )r H , (3.27) 



^> : -iGV(vp q )r H , (3.28) 

with a momentum conservation factor (27r) 4 £ 4 (At> — k — p q ), where Tp = 275 for pseudoscalar 
(H — P), T v — — for vector mesons (H = V), and p g is the momentum of the light degrees 
of freedom. 

(ii) The internal line propagators for the heavy quark and the light antiquark are, 



(3.29) 



2(v-k + ie) 



-p a + m q 



'^-r^- > (3-30) 



-p q Pq-m z q + ie 

respectively, where k is the residual momentum of the heavy quark, and m q the constituent 
mass of the light antiquark. 

(iii) For internal lines, integrate over the internal four-momenta, 

d 4 k r d A p Q 

and 



(2tt) 4 ' J (2tt) 

for heavy and light quarks, respectively. Also there is a factor of (—1) for each fermion loop. 

(iv) For all other vertices that do not attach to the bound states, the corresponding 
diagrammatic rules are standard from the the conventional field theory formulation. Most 
of these vertices mainly come from C mQ , hence the corresponding 1 / mq corrections can be 
obtained from the standard perturbation field theory (see an explicit example in the next 
section). 

These are the Feynman rules needed for subsequent calculations in this effective field 
theory. 



IV. HEAVY MESON PROPERTIES IN HEAVY QUARK LIMIT 

In this section, we evaluate the basic heavy meson properties in the heavy quark limit 
within the present framework. These include the Isgur-Wise function, the decay constants 
and the axial- vector coupling constants of heavy mesons. 
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A. Isgur-Wise function 



In the heavy quark limit, the transition matrix elements of B — > D, and B — ► D* are 
given by 



(D(v f )\fi v/ Th b v \B(v)) and (D(v', e*)\h v ,Th v \B(v)} , (4.1) 

where T is a Dirac 7-matrix from the electroweak current of heavy quarks. Using the 
Feynman rules given in Sec. IIIC, the hadronic matrix elements of B — ► D and B — > D* 
decays turn out to be (see Fig. 2(a)) 

{D(v')%Th h v \B{v)) = Tr{r P (^)r(l±^)r P .M} , (4.2) 
(D{v\e*)\X,Th b v \B{v)) =Tr{r v (^)Y(^)T P M}, (4.3) 



or 



(H\v')\h' v/ Th v \H(v)} = Tr{T H \^)r(^)r H M} , (4.4) 
for the general heavy meson decay process H — > H', where 

M _ iC 2 f ^Vi _ gV • PiM v ■ Pi) tg-m q 

J (27r) 4 (A - v ■ p q + ie) (A - v ■ p q + ie) p\ — m 2 + ie 1 

which is actually the transition matrix element of the light antiquark (brown muck). 
Since M. is fully co variant, it has the form 

M = A + Byt + C/, (4.6) 

where the coefficients A,B,C,D can be easily determined to be 

A = -iG 2 ( P^ iT - - .Z {V '' Pq)nV '. P i rr rn q , (4.7) 



B = iG 



(27r) 4 (p 2 - m 2 + ie)(A - v ■ p q + ie)(A — v'-p q + ie) 

d% V*(v'-p q )V(vp q ) 

(2tt) 4 (p2 _ m 2 _|_ j e )(X - v ■ p q + ie)(A -v'-p q + ie) 



C = iG' 



x i|(^±^ + (^r^i (48) 

2{(l + vv') (1- vv') y K J 

d% ^*(v' ■p^jvp,) 

(27r) 4 (p 2 — m 2 + ie) (A — v ■ p q + ie) (A — v' ■ p q + ie) 
U {v + v').p {v-v').p \ 

2{{l + vv') {1-v-v') j K J 

Then Eq. ( |4.4| ) can be simplified to 

(H'(v')\W v ,Th v \H(v)) = -Z( v . v ')Tr{r H {^)r(^)r H }, (4.10) 
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where £(v ■ v') is the Isgur-Wise function 

Hv ■ v') = -(A-B-C) 
d A p q 



m q + (v + v') ■ p q / (1 + v ■ v') 



(p 2 - m 2 + ie)(A - v ■ p q + ie)(A - v' ■ p q + ie) 
At zero recoil v ■ v' — 1, i.e., i>' = 



v -pq + rriq 



(p 2 - m 2 q + ie) ( A - v ■ p q + ie) 



1 



(4.11) 



(4.12) 



where we have used the result ( |3.21|) . This is the well-known normalization of the Isgur-Wise 
function at the zero-recoil point, as dictated by HQS. 



B. Decay constants in the heavy quark limit 

The decay constants of pseudoscalar and vector heavy mesons are defined by (0|A M | P(p)) 
= ifpp^ and {0\V^\V(p,e)) = fyMye^ respectively, where = q^^^^Q and = q^y^Q. 
In the heavy quark limit, the meson decay constants are redefined as 

(0\<?f'y s h v \P(v)) = iFpv» , (0\qYh v \V(v,e)) = F v e^ (4.13) 



where in the heavy quark limit, \P(p)) = y/Mp\P(v)) and likewise for \V(p,e)), therefore, 



F P = f P ^M P , F V = fy^My. (4.14) 

HQS demands that F v = F P . 

Now, using the Feynman rules of the effective theory, it is very simple to evaluate the 
above matrix elements (see Fig. 2(b)): 

(0\q^ l5 h v \P(v)) = Tr{ 7 ^7 5 ^r P ^ 1 } , (4.15) 

(0\qYh v \V(v,e)) =Tr{f^-TvMi}, (4.16) 



where 



M l = iGy/N e f^* i V(v.p 9 



dPq m , \ i> q -m. 



and 



(2vr) 4 (p 2 q - m\ + ie) (A - v ■ p q + ie) 

Ai + Btjf, (4.17) 



Ai = -iGsfN c f £%*(v ■ Pq )- — ^ — , (4.18) 

v J (27r) 4 (p 2 - m 2 + ie){A - v ■ p q + ie) 

Bl = lG JK f ^%nv-p q )- 2 ^— —t • (4.19) 

v J (27r) 4 (p 2 - m 2 + ie)(A - v ■ p q + ie) 
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Here N c = 3 is the number of colors. Thus, it is easily found: 
F P = -2{A X - B x ) 



= 2iG [Nc ( q <L'(v ■ p ) v-p q + m q 

v °J (27r) 4 q (p^ — rriq + ie)(A — v ■ p q + ie) 

= F v , (4.20) 
as expected from HQS. 



C. Axial- vector coupling constants in the heavy quark limit 

The axial- vector coupling constants, / and g, are defined by the transition matrix ele- 
ments of V — > Pn and V — ► V rr in the soft pion limit : 

(P(v)\A°\V(v,e)) , (V'(v,e'*)\A«\V(v : e)}, (4.21) 

where A a = <? 4^7^75 ? is the light quark axial- vector current. In the heavy quark limit, 

(P(v)\q-A a \V(v,e)) = it e .q, (4.22) 
iy'{v ) e'*)\q-A a \V{v ) e)) = ige^q/lv p e a , (4.23) 

where q ~ is the momentum carried by the soft-pion, and the SU(3) flavor matrix element 
X'h'^Xh has been omitted in the above expressions. HQS predicts that f — 2g 

Diagrammatically, the above two matrix elements are represented by Fig. 2c, from which 
one can straightforwardly write down these matrix elements: 



(P(v)\q ■ A\V(v, c)> = iTr{r P ^r y ^ 3 } , 

(V'(v, e'*)\q ■ A\V(v, e)> = ):Tr{T* v ^T v M 3 } , (4.24) 



or simply 



>H'(v)\q ■ A\H(v)) = It t {t h ,^±It h M 3 } , (4.25) 



where 



M3 lU J (2tt)^ [V Pg)l (p>- m l + iey(A-v.p q + 

= [A 3 v ■ q + B 3 4+ C 3 v • q j + D 3 4 y) T5 . (4.26) 

and the coefficients are given by 



13 



J (27r) 4 9 (p^ _ m 2 _|_ ^ e )2(A — v ■ Pg + ie) 

^ -n* f dA V« \*,( m2 2 4(?;-p, ? ) 2 -p 2 , 

Cs = «G / — — --^ 2 ^ . wT " -77, (4.29) 

J (,27rj 4 3 (p* — m A + ze)^(A — t> • p q + ie) 

D 3 = -A 3 . (4.30) 
Then, it is easily seen that Eq. (|4.25 ) can be simplified as 



(H'(v)\q ■ A\H{v)) = 9 -Tr{T Hl ^-T H #y 5 } , (4.31) 



2^2 

with f — 2g and 

g=-(B 3 + D 3 ) 

_.G 2 f d% hTf ^ ^ u2 (p 2 -m 2 q ) + 2(vp q + m q )(vp q + 2m q 



3 7 (2tt) 4 ' 1 ^ (p2_ m 2 + , e)2 ( A _, ; . p(? + , e ) ■ ^ ° ) 



V. THE DETERMINATION OF HEAVY MESON MASSES 

In this section, we determine the heavy meson masses up to order 1/mq within the 
effective theory. 

In the heavy quark limit, the heavy meson masses can be written as Mh = tuq + A. In 
this limit, the pseudoscalar and vector mesons are degenerate. The correction to Mm comes 
mainly from the leading HQS-breaking 1/mg corrections [see Eq. ( |3.4j) 1: 

d = h v (tD ± ) 2 h v + ^h v a^G^h v = 1 + 2 , (5.1) 

where = §[7^,7^] and G^ v = D v ]. With these I/thq corrections included, the 

heavy meson masses can be expressed as 

M H = m Q + A- -^—(X 1 + d H X 2 ) , (5.2) 

where Ai comes from 0\ and <iffA 2 from (9 2 (see Fig. 3). A, Ai and A 2 are the basic parameters 
in HQET. The parameter A is the residual mass of heavy mesons in the heavy quark limit 
and is associated with the dynamical mass of the brown muck j3|, Ai parametrizes the 
common mass shift, and A2 describes the effect of the hyperfine mass splitting. In the rest 
frame of the heavy meson, the Clebsch- Gordon coefficient du in Eq. Q5.2| ) is conventionally 
normalized in such a way that 

d H = -(H(v)\4S Q -S £ \H(v)) 

= -2[5 tot (5 tot + 1) - S Q (S Q + 1) - S e (S e + 1)], (5.3) 
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where Sq (Sg) is the spin operator of the heavy quark (light degrees of freedom). Therefore, 
d>H = —1 for vector (1~) mesons, and dn = 3 for pseudoscalar (0~) mesons. 

Using the Feynman diagrammatic rules of the effective theory, we can immediately write 
down the expressions for Ai and o?#A 2 for heavy mesons as follows: 

A 1 = Tr|r // ^r^ 2 |, (5.4) 
d H \ 2 = Trlr H ^a^^-r H M^\ , (5.5) 



with 



J (27r) 4 (jp* - m 2 , + ie) (A - v ■ p q + iey 



^ ■ (gig + P2q)v fl ~ (gig + P2g) M 

(A - u • pi ? + ie)(A - w • p 2g + ie) 
^, ^""'. 7- /""J";. ■ (5.6) 

£>2g - m q + *C Pig ~ m q + ^ 



(2tt) 4 (2tt) 4 

(Pig " P2g) M 



X -=- 



(A - v • pig + ie)(A - w • p 2g + ie) 

X D ~ o | ■ lv'~2 o | ■ ; I 5 - 7 ) 



where Cf is a color factor, C/ = ^ 1 =4/3 for iV c = 3, and D^ u is the gluon propagator 

given by D^ v = in Feynman gauge (q = p lq - p 2q ) ■ 

The mass shift parameter Ai can be further simplified to 

4^ 0^2 _(,, ~ \2 



^P, |, Tf / M2 2(p2-(V-Pg) 2 )(vPg+mg) 



(27r) 4 (p^ — m 2 + ie)(A — t> • pg + ie) 5 

-2g 2 C f G 2 f d4piq dAp2q ' ' 



(2tt) 4 (2tt) 4 (A - u • p^ + ie)(A - v ■ p 2q + ie) 
1 



x 



{(Pig - P2q) 2 + ie)(p 2 g - m 2 + ie)(p|g - m 2 + ie) 

X { (V ■ Plq + m q ) [(pig + p 2 g) • P2g - V • (pig + p 2 g)^ • P 2 g] 

+ (u • p 2q + m q ) [(pig + p 2 g) • pig - v ■ (pig + p 2( > • Pig] } • (5.8) 
The hyperfine mass splitting parameter A 2 can also be simplified. By using the identity 

^^" = 0, (5.9) 



15 



we find that only the antisymmetric component of Ai 2 contributes to A2. Thus we can let 

--I^C- (5.10) 
With the use of T H = — jfT H , it is not difficult to find that ( is given by 
t 4 % n n 2 f d A Pi q d% q • p lq )V*{v ■ p 2q ) 



r s 1 J (2ir)*(2ny(A-vp lq + te)(A-v-p 2q + te) 

1 



((P?i - P2 q ) 2 + ie)(pl q -m 2 q + ie)(p% - m\ + ie) 

X {(V ■ Pig + m q )[(Plq ~ P2q) ' P2q ~ V • (pi, - p 2f > ' P2g] 

- (V -p 2 g +m q )[(p lq -p 2q ) -Plq-V- (Pig ~ P2q)v ■ Pl q }\ . (5.11) 



Since 



,( r „l±^i±2^1 f *r ; -*. ( , 12) 



it follows that 



Therefore, 



^ ^ _ ^ x 1 3 for pseudoscalar mesons, ,g ^ 

" 1—1 for vector mesons. ' 



A 2 = C- (5.14) 

The above results are valid in any Lorentz frame. In other words, we obtain the result of 
Eq. ( |5.3| ) without setting the heavy meson in the rest frame. From Eq. (|5.2|) , we obtain the 
hyperfine mass splitting, 

2A 

AM VP = M v - M P = — . (5.15) 

m Q 

As we can see from Eq. (~ |5.8|) , Ai consists of two contributions: one is the kinetic energy 
of heavy quarks [the first term in Eq. (BT8[)] , and the other is the effect of the chromo-electric 
interaction between the heavy quark and light degrees of freedom. For convenience, we 
redefine 

Ai = -(/?) + a s Ai , A 2 = a s \ 2 , (5.16) 

where a s = g^/Ait, and (k 2 ), Ai and A 2 are denoted as the kinetic energy, the chromo-electric 
and chromo-magnetic interaction parameters, respectively. These parameters depend only 
on the structure function ty(v • p q ). 
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VI. COMPARISON OF THE EFFECTIVE THEORY WITH THE CO VARIANT 

LIGHT-FRONT MODEL 



In the previous sections we have developed an effective field theory to describe heavy 
mesons as composite particles, and systematically evaluated the heavy meson properties in 
the heavy quark limit and the basic HQET parameters. In this section, we will compare 
the effective theory with the covariant light-front model of heavy mesons that we have 
constructed recently |14}| , and show the source of difficulty for developing a fully covariant 
reformulation of the currently used light-front quark model in the literature. 



A. Covariant light-front quark model as a special case 

So far we have not specified the form of the heavy meson structure function \l/(t> • p q ). 
From the constraint that it forbids on-mass-shell decay of the heavy meson into Qq, a 
possible form is given by 

$(v-p q ) = (K-v -p q )<p(v-p q ), (6.1) 

where the function ip(v-p q ) is regular at v-p q = A. We shall further assume that (1) <p(v-p g ) 
is analytic except for isolated singularities when continued into the complex plane, and (2) 
it vanishes "fast enough" as \v ■ p q \ — > oo. These two conditions allow us to perform the 
dp e — integrations in the expressions derived earlier by Cauchy's Theorem. Thus by closing 
the contours in the lower-half complex p e — plane, we obtain 

avv')=iG 2 J 0^^W-pMvp q ) 

m q + (v + v') ■ p q /(l + v ■ v') 



(P q -m 2 + ie) 



drp q 



(27T) 



x 



rn 



+ (v + v')-p q /(l + v-v') , (6.2) 



2 / v / d % t ^v-p q + m q 



and 



(2vr) 4 ^ v rq 'pl-ml + ie 
■2G\/X. I j^(27r)5(pl-ml)e(p G q Mv-p q )(v-p q + m q ), (6.3) 



G- 2 = J ^(27r) ( 5(^-m^(p e )(,;-p (? + m g )|^-^)| 2 . (6.4) 



Now if we take 



y/v-p q + m q 
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where &lf{v ■ p q ) is a normalized light-front wave function, and let p e = p , then these 



results are exactly the same as those obtained in the covariant light-front quark model [14 



Thus we see that, for these quantities, the covariant light-front quark model results can be 
recovered here with a special choice of tp{y ■ p q ). 



B. Beyond the covariant light-front quark model 

Next, we shall consider the mass parameters Ai and A2 (see Sec. V). As before, using 
Eq. ( |6.1|) and after performing the p e -integrals, we obtain 

Ax = 2G 2 / ^i(^)5(pl - ml)6(pl)\rtv ■ p q )\ 2 [m 2 q - (v ■ p q f] (v ■ Pq + m q ) 
f (2n)5(pl q - m 2 q )e(pt q )(27r)5(pl q - m 2 q )9{p% q ) 



x 



1 (Pig - P2g) 2 + ^ 

(2n)6(p 2 lq - my(pi q )(2n)6((p lq -p 2q ) 2 )9(p%-pt q 



+ 



Vl q ~ ™ q 

(27r)S(p 2 2q - m 2 q )e(p%)(2n)5((p lq -p 2q ) 2 )6(pi q -p% 



Vi q ~m 2 j 

X { (V ■ Plq + m q ) [(pi q + p 2q ) -P2q-V- (p lq + p 2q )V • p 2q ] 

+ (v ■ p 2q + m q ) [(p lq + p 2q ) -piq-v- (p lq + p 2q )v ■ p lq ] } , (6.6) 



x 



+ 



(2vr) 4 (2tt) 

f (2n)5(p 2 lq - m 2 q )9(pt q )(2n)5(p 2 q - m 2 q )9(p%) 

\ (Plq ~ P2qf + It 

(27r)5(p 2 q - m 2 q )6(pt q )(2n)5(( Plq - p 2q ?)6{pl q - p\ q 



Vl q ~ ™ q 

(27T)5(p 2 q - m 2 q )6(p%)(2n)5((p lq -p 2q ) 2 )6(pi q -pl q 



Viq ~m 2 q j 
x { (v ■ Pig + m q ) [{pi q - p 2q ) -p 2q -v- (p lq - p 2q )v ■ p 2q ] 

- (v ■ (p 2q + m q ) [(p lq - p 2q ) -pig-v- (p lq - p 2q )v ■ p lq ] } . (6.7) 

In each of the above expressions, the first delta-function term inside the bracket comes 
from on-mass-shell light antiquarks. With (p(v ■ p q ) given by Eq. (|6.5|) , one can easily check 
that this contribution is exactly what one would obtain in the covariant light-front quark 
model [Q (where the expression for Ai is not explicitly displayed). The rest of the terms 
are due to the off-mass-shell contributions of the antiquarks, they correspond to the so- 
called Z-diagram contributions which cannot be obtained in any conventional quark model 
formulations. 
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As another example, we examine the axial-vector coupling constant g, which involves a 
purely light-quark current. To obtain a covariant result, one has to keep all the time-ordered 
(including the so-called Z-diagram) contributions, which is obviously beyond the scope of the 
light-front quark model |Tl| . However this is not a problem in the present framework, since 
Feynman diagrams naturally contain all possible time-orderings. As a result, covariance is 
automatically preserved in deriving all the results in the previous sections. 

Now substituting Eq. QS.1| ) into Eq. ( |4.32[ ), and integrating over p~ , we obtain 



9 = Y 



G 2 r dp+dp q± ( 



(27r) 3 2p+ \ 



¥(V-Pq)[ 



,(vp q + 2m q )(vp q + m q ) (A - v ■ p q ) 



X 



x 



2v ■ p q + 3m q + X + (v ■ p q + 2m q ) (v ■ p q + m q ) 



d 



d(v-p q )_ 



X 



\<P(V-Pq)\' 



(6.8) 



where X = p^ /v + and p 2 q 



m 



9" 



Obviously, this expression cannot be obtained from any 



type of light-front quark models because it involves a derivative of the wave function, which 
is very unusual and cannot be simply related to some matrix elements of hadronic bound 
states in quark models. However it is interesting to note that the first term in Eq. ( |6.8| ) 
is what one would get for g if we naively calculate it in the light-front quark model ||14|| . 
Therefore the second term in ( |6.8[ ) should be the Z-diagram contribution not present in the 
light-front quark model formulation. 

For completeness, we also display the resulting expression for g by integrating over p° in 



Eq. (FQ21) : 



G 2 
9 = T 



dp 3 ( 2 (vpq + 2m q )(v -pg + mq) (A-v-pq) 

wr ( " ' Pq)l — 



X 



+X + (v -p q + 2m q )(vpq + m,)( - \ + 




where X = p° jv and p 2 



mz 



P° 8(V ■ p q ) 



Thus we have shown that, by the specific choice of Eq. (|6.5|) , our effective field theory 
can reproduces the covariant light-front quark model |l4j] results for those quantities (£(v • 
v'),Fp, and Fy) which do not have Z-diagram contributions. For other quantities (Ai,A2, 
and g), the extra Z-diagram contributions missed in the light-front quark model can be 
explicitly identified here. Hence the effective field theory presented here has the simplicity 
of a conventional quark model, and is fully covariant at the same time. 



C. Choice of ^(v ■ p q ) 

Now an important question is what kind of ^(v-p q ) will allow a proper analytic continua- 
tion into the complex plane, so that the p e -integrals can be evaluated by Cauchy's Theorem, 
and comparisons can be made with results obtained in the quark model. In the light-front 
formulation, there are three type of mesonic wave functions which have been widely used in 
the literature. One of them, the so-called BS W wave function (or BSW model) has already 
been ruled out since it explicitly breaks the relativistic covariance and thereby leads to in- 
consistent results in HQET, as we have shown in recent publications • The other two 
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wave functions, namely, the Gaussian-type wave function and the invariant light-front mass 
wave function, have the following forms in the heavy quark limit [14]: 



^f(v ■ p q ) = AT Gv ^7^exp I - [( v . Pq f - m 2] I ? ( 6 .io) 
< F (v-p q )=Af M expl- V -^\. (6.11) 



These two wave functions apparently have a covariant form [TJj]. However, in the light-front 



p 1 -(-ttt, 2 ^ 

quark model formulation, the light antiquark is always on the mass shell: p~ = "* x + "' . It 



Vq 



is not obvious how to extend these wave functions to be used in a 4-dimensional covariant 
framework. If p q were allowed to go off-shell in these wave functions, then $>ip(v ■ p q ) 
would be unbounded when \p e \ — > oo, whereas $1f( v ' Pq) w °uld also be unbounded when 
analytically continued into to the complex plane. To fix these problems, one could instead 
use \v ■ p q \ in the above expressions, but then both wave functions would not be analytic in 
the complex plane. In any case, we find that it is not possible to modify the exponential 
form so that the two conditions stated below Eq. fl6.1|) are satisfied. 

Besides the exponential-type functions, we can also consider the so-called Lorentzian- 
type dependence for ^{y ■ p q ). In the heavy quark limit, we may write ^f(v ■ p q ) as 

* n (v-p q )= A ~ V ' Pq .. n , (6-12) 
[v -p q + u — ie) n 

where an — ie has been inserted in the denominator so that *f> n (v-p q ) is analytic in the lower 
half complex-^ 6 plane. Consequently by closing the integration contours in the lower half 
complex-^ 6 plane, we will not pick up contributions from the poles of ^ n (v ■ p q ). It is easy 
to check that this choice of *f?(v ■ p q ) corresponds to a covariant light-front wave function of 
the form 

n F {v-p q )=M L ^- p ^ m l n>2, (6.13) 

[V ■ p q + UJ) n 

where Ml is a normalization constant, and p 2 = ml. Here we require that n > 2 to ensure 
the finiteness of all the integrals. 

In the following, we shall present some numerical analyses to further examine the pre- 
dictive power of the effective theory with the choice Eq. ( |6.12[ ) for ty(v ■ p q ). 



VII. NUMERICAL CALCULATION AND DISCUSSIONS 

A. Short-distance corrections 

Before performing the numerical computation, it should be stressed that all the quantities 
we have evaluated in the effective field theory are mainly governed by the long-distance 
physics of HQET. That is, a renormalization scale \i (Aq C d « A* << m o) has been 
implicitly employed in the structure function ty(v ■ p q ). It is necessary to take into account 
short- distance QCD corrections to match with the full QCD description: 
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((9> Q CD = C O (/l)(0o(/i)>HQET + 



CM 



HQET 



+ 



(7.1) 



where the Wilson coefficients Cj(/i), which account for the short- distance corrections, have 
been computed in the literature ||. 

Explicitly, the B— meson decay constant J'b is given by 



h 



a s (m b ) 



-i -6/(33-2Af / ) 



Fp{h) 



(7.2) 



where Fp(fi) is the decay constant defined in Sec. IVB, and Nf = 4 is the number of quark 
flavors with mass less than m&. 

The Isgur-Wise function discussed in Sec. IV. A is also defined at the renormalization 
scale fi. Its relation with the //-independent physical observable, for example, the form 
factor f+{q 2 ) [or F^q 2 )} is given by |21| 



U(Q 2 



M B + M D 



2^M B M D 



a s (m c ) 



-i -6/25 



-i a L (v-v') 



(7.3) 



where 



a L (u) 
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[ur(u) - 1], 



vV^i 



ln(u; + Vuj 2 - 1). 



(7.4) 



The axial-vector coupling constant g does not receive short-distance corrections since it 
involves a light-quark current which is partially conserved. 

Although the physical heavy meson mass is scale independent, individual contributions 
to it are modified by short-distance corrections. Thus Eq. (|5.2j ), written separately for 
pseudoscalar and vector mesons, becomes 



m Q (/i) + A(/i) 
m Q (fi) + A(fi) 



1 



1 



Ai(/i) + 3a 2 (/i)A 2 (/i) 
Ai(/i) - a 2 (/x)A 2 (/x) 



(7.5) 
(7.6) 



where 



a 2 (/x) 



a s (/f) 



-i -9/(33-2^) 



(7.7) 



Note that all the basic nonperturbative parameters in HQET, i.e. m^/i), A(/a), Ai(/x) and 
A 2 (/a), are //-dependent; only the physical masses are independent of the renormalization 
scale \i. 

It should be stressed that in the above expressions, we have distinguished the heavy 
quark mass in the I/tuq expansion which is the running heavy quark mass defined at the 
scale /i, and the mass scale in the strong coupling constant Q s (mg) set at the heavy quark 
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pole mass: a s (mg le ). In the literature, there are many discussions on the ambiguities in the 
definition of the heavy quark mass due to the presence of the so-called renormalons [J22[ |. We 
believe that the above distinction of heavy quark masses in our formulation is theoretically 
consistent. 

Using Eq.( |5.16|) , Eqs. ( |7.5| - |7^ ) can be rewritten as 



M H = mq(ji) + A(/i) + - — ; {(k 2 ) - a a (n) \Xi + 3a 2 (^)A 2 l ) , 
M H * = m Q (fi) + A(n) + - — f (k 2 ) - a s (//)[Ai - a 2 (//)A 2 l) . 



(7.8) 
(7.9) 



Note that (/c 2 ), Ai and A 2 depend only on the structure function ty(v-p q ), and therefore also 
on the scale fi. The heavy meson mass splitting between pseudoscalar and vector mesons is 
simply given by 



M H *-M H = 2 a ^ )a / M \ 2 . 



(7.10) 



B. Numerical analysis 

For the sake of demonstration, we first set \i = m^ ole ~ 4.89 GeV (the pole mass) 
27 1, and see if a consistent set of parameters can be found. At this scale, short- distance 
corrections are not present. The light quark mass at this scale is given by the current quark 
mass which is about a few MeV for up and down quarks. We take m Ut d(ml° lc ) ~ 5 MeV. 
Then the parameter u in the structure functions can be fixed by the decay constant via 
Eq. (pL20p . Taking 



f B = 0.180 GeV, (7.11) 

the values of lo for different structure functions ^f n (v ■ p q ) are listed in Table I. Since the 
strong coupling constant a s is also known at fi = m^ ole , knowing u>, we can predict various 
heavy meson properties obtained in the effective theory at this scale. 

We first calculate the vector-pseudoscalar B meson mass difference which is determined 
by A 2 . As /i = ml° lc ^> A QCD , a s (m^ ole ) can be determined from a s (M z ) by the perturbative 
evolution equation (at the one-loop level in the MS scheme), 

a (m p ° lc ) = Qs(Mz) ~ 22 (7 12) 

b l + a s (M z )(3 \n[(m p b ole ) 2 /M 2 ]/(An) ' 

where f3 = 11 - |JV> = 11 - 8/3 for N f = 4, M z = 91.19 GeV, and a s (M z ) = 0.119 from 
experimental fits. Thus, the B* — B mass splitting at the scale \x = ml° le in our calculation 
is given by 



a. Am 



pole\ 
b 



AM b *b = M B * -M B = 2 _) " ' A 2 , (7.13) 

m h [ml ) 
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with the MS 6-quark mass at this scale: m&(m£ ole ) = 4.339 GeV p7j- The result is also 
listed in Table I, from which it is evident our results do not fit the experimental value of 
AMb*b — 0.046 GeV ||23|| . This numerical analysis indicates that the real scale \x implicitly 
used in the effective theory should be lower than mj° le . 



Table I. The parameter to in the structure functions \& 
fs = 180 GeV at the scale /j, = m£ ole 
of AMb*b and (k 2 ) becomes insignificant 



■p q ) m 



in; ' and the resulting AMb*b and (fc 



T2] )] with different n values fitted to 
; ). With n larger than 12, the change 



n 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


w (GeV) 


0.100 


0.374 


0.689 


1.011 


1.335 


1.660 


1.986 


2.312 


2.638 


2.965 


AM B *B (GeV) 


0.011 


0.023 


0.028 


0.031 


0.032 


0.033 


0.034 


0.034 


0.035 


0.035 


(fc 2 ) (GeV 2 ) 


0.111 


0.274 


0.374 


0.432 


0.470 


0.496 


0.515 


0.529 


0.541 


0.550 



Of course, as is well known, the HQET works well only at the scale A.qcd << /i << m.&. 
The results shown in Table I indicate that the implicit scale fi of our effective theory is 
not close to m&. On the other hand, for perturbation expansion to work, p, should not be 
too close Aqcd either. Hence, we must choose a value of p such that the corresponding 
(XsKqcdI^q is small enough to ensure the validity of the HQET. In Table I, we have also 
calculated (k 2 ), the mean square momentum carried by the quarks in the meson. We find 
that {k 2 ) increases with increasing n, which indicates that in fact the scale p also implicitly 
increases with n. As a result, the corresponding a s would be larger for smaller n. Thus, 
in order that higher order I/tuq corrections can be ignored, we should choose a structure 
function ^ n (v ■ p q ) with a suitably larger n. 

In the following numerical analysis, we shall consider the scale with 

a s (p) < 0.5 and m q (fi) = 0.22 ~ 0.25 GeV (7.14) 

and try to find a consistent set of parameters for each ^ n (v -p q ). The condition a s {p) < 0.5 
guarantees that a s (p)AQCD/mb{p) cannot be too large. Also, at this scale, we hope that 
the short-distance corrections given in Sec. VILA are still applicable. The range of the light 
quark mass we have assumed is close to the constituent quark mass often used in relativistic 
quark models. Now we proceed as follows. For each *& n (v -p q ), to is fixed by fs = 180 GeV 
through Eq. ( |7.2| ) . We can then predict the HQET parameters Ai and A 2 . Subsequently, 
rrib{n), A(ju), and a s (fi) are determined by fitting the experimental AMb*b — 0.046 GeV 
f2| and the experimental E-meson mass M B = 5.279 GeV |23| via Eq. ( |7TTU|) and Eq. (fT8|), 
respectively. We found that for n < 3, ^ n (v ■ p q ) allows no consistent fit to the experimental 
data. Only with n > 4, can we have good descriptions for various B meson properties with 
small I/nib corrections. The results are summarized in Table II. 



Table II. The parameter lu in the structure functions ^ n (v ■ p q ) [( 6.1 2| ) ] and the bottom quark mass 



mb fitted to fs — 180 GeV and AMb*b — 0.046 GeV with a suitable giving scale by a s (fj,) and m q (n) for 
each giving n in Eq. ( |6.12| ) , and all the predicted nonperturbative HQET parameters in the effective theory, 
where m qj uj, nib and A are in unit GeV, and (A: 2 ), a s Ai, Ai and A2 in unit GeV 2 . 
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*n(«*Pff) 




m q (p) 






<fc 2 > 


a s Xi 


-A x 


A 2 


A 


m c {p) 


9 


n = 5 


0.457 


0.251 


0.405 


4.907 


0.292 


0.254 


0.038 


0.087 


0.333 


1.626 


0.380 


n = 6 


0.431 


0.248 


0.671 


4.901 


0.329 


0.265 


0.064 


0.089 


0.337 


1.614 


0.329 


n = 7 


0.417 


0.245 


0.937 


4.894 


0.352 


0.271 


0.081 


0.090 


0.342 


1.604 


0.328 


n = 8 


0.408 


0.242 


1.205 


4.888 


0.368 


0.275 


0.093 


0.090 


0.347 


1.595 


0.337 


n = 9 


0.401 


0.238 


1.477 


4.880 


0.381 


0.278 


0.103 


0.091 


0.354 


1.584 


0.372 


n = 10 


0.395 


0.234 


1.752 


4.869 


0.391 


0.280 


0.111 


0.091 


0.364 


1.572 


0.428 


n = 11 


0.390 


0.230 


2.029 


4.859 


0.399 


0.282 


0.117 


0.091 


0.373 


1.560 


0.479 


71 = 12 


0.385 


0.226 


2.310 


4.844 


0.406 


0.283 


0.123 


0.091 


0.389 


1.542 


0.556 


exp. 








4.89** 








0.090* 


0.33 - 0.45 Hi 


1.59** 


< 0.7 Hi 



*In the literature, the usual experimental value of A2 — 0.12 GeV 2 is obtained by taking nib = Mb approx- 
imately. Here the experimental A2 — 0.090 GeV 2 is obtained by using ml° lc = 4.89 GeV, the 6-quark pole 
mass, in order to make a consistent comparison with the theoretical calculation used by Eq. ( 5.1 5| ) with the 
short-distance correction. 

** There are no experimental data for the constituent quark masses mj, and m c . Here we list the bottom 
and charm quark pole masses p7[ 

From Table II, we see that (k 2 ) increases with increasing n. As noted earlier, this 
indicates that the scale \i in the phenomenological structure function ty n (v -p q ) also goes up 
with n. Then to be consistent, the running strong coupling constant a s {p) and the running 
quark masses m q (n) and mb(//) should decrease with increasing n. This expected behavior is 
clearly shown in Table II. From table II, it is also interesting to see that the resulting 6-quark 
mass is in between the constituent mass (~ 4.8 GeV) used in various model calculations 
and the pole mass (4.89 ± 0.05 GeV). Furthermore, we also find that the hyperfine mass 
splitting parameter A 2 is quite stable and insensitive to the function ty(v ■ p q ). This is not 
unexpected since A2 is directly related to the heavy meson mass splitting which is well 
measured. The mass shift parameter Ai is very small in our calculation. In this work, Ai 
receives two contributions: one is the heavy quark kinetic energy ({k 2 ) ~ 0.3 ~ 0.4 GeV 2 for 
n > 5) which is a negative contribution to Ai [see Eq. (|5.16 )1, and the other is the chro mo- 



electric interaction between the heavy quark and light degrees of freedom characterized by 
the parameter Ai which is close to 0.3 GeV 2 and is a positive contribution. Therefore, the 
kinetic energy is largely compensated by the chromo-electric interaction energy, leading to 
a small and negative X\: Ai = 0.04 ~ 0.11 GeV 2 for n = 5 ~ 10. Then a fit to the physical 
masses of the B mesons yield the nonperturbative HQET parameter A: A ~ 0.33 ~ 0.36 
GeV, which is consistent with results from HQET analyses of the semi-inclusive B decay 
data ||26|1 . These numerical results provide a vote of confidence to the reliability of the 
effective theory proposed in this work. 

For a further consistency check with HQET and HQS, we use the above results to cal- 
culate the D meson observables. The D meson mass is given by 

M fl = mc (/i)+A(/ i )-A^-^AM D . D , (7.15) 

where A and A12 are the same for both the B and D mesons. Then using the experimental 
D mass data [p3] 
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M D = 1.8645 GeV, AM D » D = M D * - M D = 0.142 GeV, 



(7.16) 



we can uniquely determine the charmed quark mass. We get (see Table II) 

m c (/i) = 1.57 ~ 1.62 GeV (n = 5 ~ 10), (7.17) 

which is close to the m c used in various quark model calculations, and also the pole mass 
mP° le = 1.59 ± 0.02 GeV |37J . 



Next we proceed to calculate the Isgur-Wise function £(v ■ v') [cf. Eq. ( |6.2| )1, its slope 
parameter at the zero-recoil point 

P 2 = -£'(1), (7-18) 



and the axial- vector coupling constant g from Eq. ( |6.8|) or Eq. (|6.9| ). The Isgur-Wise 



function depends on the parameter tu and the light quark mass m q , and the axial- vector 
coupling constant also depends on A. Hence, once u and A are determined for a given m q , 
£(t> ■ v'), p 2 , and g can all be predicted. In Fig. 4, we plot the Isgur-Wise function as a 



function of v ■ v'. It is remarkable to see that the structure functions ^f n (v ■ p q ) (|6.12| ) with 
different n give very similar results. In particular, the Isgur-Wise function obtained from 
v|/ n=10 is almost identical to that form the light-front wave functions (|6.10|) and ( |6.11| ). The 



slope parameter of the Isgur-Wise function at the zero-recoil point is found to be: 

p 2 = l.l~1.2 (n = 5~10), (7.19) 

which is also consistent with other theoretical analyses ||. 

Having determined A from the heavy meson mass calculations, we can now predict the 
axial-vector coupling constant g. As listed in Table II, we have 

g = 0.33 ~ 0.43 (n = 5 ~ 10) , (7.20) 



which is consistent with the experimental constraint g < 0.7 0,^] derived from D* — > D+ty 
decays. Our result is also close to the QCD sum rule results, which tend to concentrate within 
the range 0.2 < g < 0.4. 

VIII. CONCLUSIONS AND PERSPECTIVE 

In this paper, we have given a detail account of a recently proposed field theoretical 
description of heavy mesons |TJ] . The effective theory incorporates heavy quark symmetry 



and the heavy quark effective theory, from which a natural realization of heavy mesons in the 
heavy quark limit as a composite particle of the reduced heavy quark coupled with a brown 
muck of light degrees of freedom is provided. This theory is fully covariant, so that Feynman 
diagrammatic techniques can be use to carry out perturbative calculations. Moreover the 
effective theory preserves the simplicity of a conventional quark model, in fact, for those 
quantities which do not have the so-called Z— diagram contributions, light-front quark model 
results can be reproduced as a special case in the present theory. Thus this theory provides 
a link between the fundamental QCD and the phenomenologically successful quark model, 
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so that the difficult subject of hadronic bound state physics can be quantitatively studied 
in a covariant framework. 

The effective theory provides a quasi-first-principles description of the heavy meson dy- 
namics. Although at present the description of the heavy mesons structure function ty(v-p q ) 
is still phenomenological in nature, it offers a systematic approach to evaluate the 1/toq cor- 
rections based on the first-principles I/tuq expansion of QCD. This resembles very much 
the situation of the QCD analysis of deep inelastic scatterings, in which the low energy dy- 
namics (described by parton distribution functions) is determined phenomenologically and 
perturbative corrections are given in a fully first-principles way. Here the situation may 
even be better since the phenomenological part is constrained by HQS and HQET, and the 
nonperturbative QCD dynamics should be much simpler in the heavy quark limit. 

As we have seen, the introduction of the structure function \I/(t> • p q ) is crucial for a 
field theoretic realization of the heavy mesons as composite particles. In fact, this structure 
function is related to the covariant wave functions of heavy meson bound states, and it 
essentially describes the brown muck structure of the light degrees of freedom inside the 
heavy mesons in the heavy quark limit. Our results show that it is not possible to extend the 
popular Gaussian-type wave functions to be used in a fully covariant formulation. However, 
we find that a Lorentzian-type function can be readily adopted in a covariant formulation 
for heavy meson structures, and numerically it gives a very good description of the physical 
properties of heavy mesons. The predictions for all the HQET parameters are consistent 
with experiment. A is about 0.33 ~ 0.39 GeV [for ^ n (v ■ p q ) with n = 5 ~ 12] which is 
consistent with results from HQET analyses of the experimental data | 26fl . The heavy quark 



masses and m c , determined by fitting to the B* — B and D* — D mass differences, are 
about 4.84 ~ 4.91 GeV and 1.54 ~ 1.63 GeV respectively, which agree with those used 
in various relativistic quark model calculations, and also their respective pole masses. We 
find that the mass shift parameter Ai in HQET is negative and very small (about -0.04 
~ —0.12) because of a large cancellation between the heavy quark kinetic energy and the 
chromo-electric interaction between the heavy quark and the light degrees of freedom. The 
Isgur-Wise function we obtained is also consistent with other calculations. Thus, we have 
self-consistently determined and predicted all the HQET parameters within the effective 
theory. 

We expect that the effective theory presented in this paper can now be applied to describe 
various heavy meson processes, such as the inclusive B decays and various exclusive B decays 
without involving light mesons. For processes involving light mesons, we must first extend 
our framework so that light mesons are also included. Finally, to solve the structure function 
\l/(f -p q ) from the fundamental theory in the heavy quark limit is an important and interesting 
challenge. We shall leave these topics to future investigations. 
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Figure Captions 



Fig. 1 Heavy-light quark scattering in (a) quark-quark coupling picture, and (b) meson- 
quark coupling picture, which determine the composite particle structure of heavy 
mesons. 

Fig. 2 Feynman diagrams for (a) the Isgur-Wise function, (b) decay constant, and (c) 
strong axial-coupling constant of heavy mesons. 

Fig. 3 Feynman diagrams for I/tuq corrections to heavy meson masses. 

Fig. 4 The Isgur-Wise functions as a function of v ■ v' that obtained from the wave 
functions *$? n (v ■ p g ), (n = 6, 8, 10) [( |6.12|) 1 and compare with the Isgur-Wise functions 



obtained from the light-front wave functions \P [( 6.10 )] and \P [ (|6.11 )] [13 . 
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FIG. 1. Heavy-light quark scattering in (a) quark-quark coupling picture, and (b) meson-quark 
coupling picture, which determine the composite particle structure of heavy mesons. 
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FIG. 2. Feynman diagrams for (a) the Isgur-Wise function, (b) decay constant, and (c) strong 
axial-coupling constant of heavy mesons. 
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Ai 



\v\ ~ (v ■ p q ) 2 } 
Av Av 




-g s T a [v ■ (pi +p 2 K - (P1+P2)"] 

--'-<T>-- 

ig s T a Y 



Av 




ig s T a -f u . 



FIG. 3. Feynman diagrams for 1/rriQ corrections to heavy meson masses. 
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FIG. 4. The Isgur-Wise functions as a function of v ■ v' that obtained from the wave functions 
^n{v ■ Pq), (n = 6,8,10) [ Q6.12D ] and compare with the Isgur-Wise functions obtained from the 
light-front wave functions * G [ (gjjj) ] and V M [ (loUl"!) ] @. 
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